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Abstract .  Constitutive  equations  are  discussed  for  a  mixture  of  any 
number  of  materials  with  elastic  and  viscous  properties  in  which  the 
constituents  may  have  different  temperatures . 

1.  Introduction 

The  recent  literature  on  the  continuum  theory  of  mixtures  contains 
rather  diverse  points  of  view  in  the  primitive  concepts  and  the  forms 
of  the  basic  field  equations.  The  basic  equations  representing  balances 
of  masses,  linear  momenta,  moments  of  momenta,  and  energy  for  each  con¬ 
stituent  together  with  equations  for  the  mixture  as  a  whole  obtained  by 
summation,  as  developed  by  Truesdell  [l,2],  Truesdell  and  Toupin  [3]  and 
a  number  of  other  writers  are  closely  related  and  are  based  on  similar 
points  of  view.  On  the  other  hand,  the  development  of  the  theory  by 
Green  and  Naghdi  [4]  and  Mills  [5]  is  based  on  different  primitive  con¬ 
cepts  and  consequently  some  of  their  resulting  equations  have  different 
forms  from  those  in  [1,3] .  However,  in  a  recent  paper  confined  to  mix¬ 
tures  with  a  single  temperature,  the  differences  and  the  relationship 
between  the  two  approaches  have  been  clarified  by  Green  and  Naghdi  1 6 J : 
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The  basic  equations  in  the  two  developments  are  equivalent,  although 

this  is  not  apparent  at  first  flight  in  the  case  of  some  of  the  equa- 
* 

tions  --  the  differences  lie  chiefly  in  the  primitive  concepts,  the 
premise  under  which  the  twc  forms  of  the  theory  are  developed,  in 
addition  to  the  interpretations  associated  with  some  of  the  quantities 
which  occur  in  the  equations.  For  further  detail  regarding  the  relation¬ 
ship  between  the  two  forms  of  the  theory,  as  well  as  additional  references 
on  the  subject,  we  refer  the  reader  to  [6]. 

In  the  present  paper,  we  are  concerned  with  a  mixture  of  any  number 
of  constituents  each  of  which  is  at  a  different  temperature.  Since  we 
now  regard  the  previous  paper  of  Green  and  Naghdi  [7]  involving  many 
temperatures  as  too  restrictive,  we  first  reconsider  here  the  forms  of 
the  basic  equations,  as  well  as  the  entropy  inequality  for  the  mixture 
as  a  whole.  Although  the  discussion  of  the  basic  equations  in  sections 
3-5  involves  some  repetition  of  known  results,  it  is  included  here  in 
order  to  have  the  basic  equations  in  line  with  our  point  of  view  [6] 
and  in  terms  of  variables  which  we  believe  provide  simple  interpretations 
when  discussing  constitutive  equations . 

In  section  6,  we  first  briefly  outline  two  procedures  for  the  dis¬ 
cussion  of  constitutive  equations  and  adopt  the  simpler  of  the  two  which 
is  slightly  more  restrictive  than  the  other.  We  then  consider  a  mixture 
of  v  constituents  (with  different  temperatures)  which  are  not  necessarily 
gases  and  which  have  viscous  and  elastic  properties.  In  particular,  we 

In  [6],  Green  and  Naghdi  have  given  a  different  interpretation  to  the 
internal  energy  in  the  energy  equation  used  by  them  in  [4]  and  have 
recast  this  equation  in  a  slightly  different  form. 
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obtain  explicit  results  for  the  "equilibrium"  values  of  various  quantities 
and  also  include  a  discussion  of  the  results  when  all  constituents  have 
the  same  temperature. 
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2 .  Notation 


We  summarize  relevant  notation  for  kinematic  quantities  and  density. 

We  consider  a  mixture  of  v  interacting  constituents,  each  of  which  is 
regarded  as  a  continuum;  we  refer  to  the  constituent  as  the  continuum 
s  (or=  1,2, . .  ,v)  .  We  assume  that  each  point  within  the  mixture  is 
occupied  simultaneously  by  all  v  constituents ,  which  are  in  motion 

relative  to  a  fixed  system  of  rectangular  Cartesian  axes.  The  position 

a, 

of  a  typical  particle  of  s  at  time  t  is  denoted  by  x^(t)  .  where 

at,  .  at,  a  at  at  .  ,  ,  , 

xJt)  =  x.(X1,  X2,  X3,  t)  ( <  t  —  t )  ,  (2.1) 

a 

and  X^  is  a  reference  position  of  the  particles  of  s^  .  All  Latin 
indices  are  tensor  indices  and  take  the  values  1,2,3,  and  the  usual 
summation  convention  is  employed.  The  Greek  letter  a  ,  as  in  (2.1),  is 
reserved  for  reference  to  the  constituent  of  the  mixture  only.  Unless 
indicated  otherwise  the  summation  convention  does  not  apply  to  Greek 
indices.  We  use  the  notation 


at 

x.  = 
i 


x“(t) 


(2.2) 


assume  that  the  particles  of  s^  (or  =  1,2, . . ,v)  all  occupy  the  same 
position  at  time  t  so  that 


1 


x. 

i 


(2.3) 


and  we  refer  to  this  position  at  time  t  as  x^  . 
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and  put 


a 

u. 

1 


(2.10) 


so  that 


I  p  u%  0  .  (2.11) 

Of=l  “  1 


We  also  define  the  operator 


(*)  .  2U  -  2L1+  v  iU 

'  '  Dt  3t  k 


(2.12) 


» 

and  observe  that 


a 

D_ 

Dt 


_D_ 

Dt 


v 

+  Z 
P=1 


fft 

P 


(va.vP)  _a.  =  i. 
'  k  V  3^  Dt 


a  3 
\  3^ 


(2.13) 


I 


6. 
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3-  Mass  and  momentum 


Using  a  fixed  surface  A  enclosing  an  arbitrary  volume  V  the 


,  + 


equation  of  mass  balance  for  s  is 

a 


it !  pa  dv  +  S  pa\  \  I  madv  ’ 


where  is  density  of  mass  production  arising  from  all  the  other 

constituents  and  n,  is  the  unit  outward  normal  to  A  .  Hence 
k 

a 

_  3pa  3  ,  or,  D  pa  jx 

ma  ~  at  +  aXj^  Pck  Vk^  “  Dt  +  Pa  dkk  ’ 


and  w-  add  the  condition 


£  m  =  0 

0=1  “ 


Next  we  consider  the  equation  of  linear  momentum  for  s  in  the 

o 


E  orm 


Alternatively  we  could  use  material  volume  elements  of  each  constituent 
which  coincide  at  time  t  . 


(3.1) 


(3.2) 


(3.3) 
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s  r  Q!  r  a  ct  r  fat 

T-  p  V.  dV  -t-  p  V.  V,  n,  aA  -  m  v.  iV 
at  J  a  1  J^aikk  Jai 

V  A  V 

f  ,  a  a*  T  of 

=  J  (Pa  F,  +  p.)dV  ►  J  nk  aki  aA  , 


a  a 

where  F.  is  external  body  force  applied  to  s  and  o,  .  is  the  stress 
1  a  hi 

tensor^+assoeiated  with  the  ingredient  s  so  that  the  total  force 

a 

a 

acting  on  s  across  A  per  unit  area  is  n,  a.  .  ■  The  term  involving 

0!  K  Kl 

v^  is  similar  to  that  used  by  Ingram  and  Eringen  [8]  and  allows  for  the 

fact  that  momentum  is  supplied  to  s  not  from  rest  but  from  all  the 

a 

other  constituents  of  the  mixture  with  some  mean  velocity  v.  which 

1 

will  depend  on  all  the  velocities  of  the  remaining  constituents,  but  it 
is  not  necessary  to  make  any  special  assumption  about  this.  The  vector 
p**  may  now  be  called  a  uiffusive  force  acting  on  s^  .  In  addition  to 
(3-4)  we  impose  the  condition 


I  (p“  ♦  m  v“)  =  0 
ar=l  1  “  1 


In  point  form  equation  (3.4)  becomes 


a  .  a  x*.  a  a  Aqt 

o,  .  ,  +■  p  (F.  -  f .  )  =  m  v.-p.-m  v 
ki.k  i  x'  er  i  i  a  i 


and  we  assume  that  p^  is  a  vector  unaltered  by  superposed  rigid  body 

motions  of  the  whoie  mixture,  apart  from  orientation.  For  convenience 
a 

we  say  that  pi  is  objective  but  other  writers  give  this  word  a  slightly 


^Truesdell  and  Toupin  [ 3l  use  a  different  order  for  the  indices  on  the 
stress  tensor  which,  in  general,  is  not  symmetric. 
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is  the  alternating  tensor  and 


3  f  a  (*  a  a  |*  Aa 

Lij  ■  5t  Jv  '  vi  xj  JA  »»  vi\xj\dA-  Jv  *«  "i  xj  av 

-  I..(s«  Fr  [,i>xj dv  •  J,  \  "k;  xj  "  -  J„  kji dv  •  (3-i2) 


In  (3.12)  X..  is  a  skew  symmetric  tensor  which  may  be  called  the  dif- 

t) 

fusive  couple  acting  on  s^  .  With  the  help  of  (3-6),  equation  (3.11) 
yields 


a 

CTCki] 


(3.13) 


Of  a 

where  Is  the  skew  symmetric  part  of  •  We  add  the  condition 


Z  Xki  “  ° 
cr-1 


(3.14) 


so  that 


CT[ki]  =  0 


(3.15) 
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r.nergy  balance  and  entropy  production  for  mixture 


Before  discussing  energy  balance  equations  for  each  constituent 
s  we  obtain  the  energy  equation  for  the  mixture.  We  assume  that 

■—  f  l(p  U  +  ^  p  v<y  v°f)dV 
at  J  V  0=1°  °  2  or  i  i' 

+  r  s(p  vvu  +  ^ p  vvv<> dA 
JAo=l  “  k  a  2  a  k  i  1  * 

■  J '<■* * s >? ’?)«*! ( s  < n-h>“  • 

V  Qf=l  A  Qf=l 


(4.1) 


.■.'here  we  take 


a  a 

t.  =  n,  <3,  . 
l  k  ki 


nk  qk 


\  =  S 


Qf=l 


V 

pr  =  Z 
Of=l 


P  r 
or  a 


(4.2) 


q  being  the  heat  flux  vector  for  s  and  q  the  heat  flux  vector 

for  the  mixture,  r^  is  heat  supply  per  unit  mass  of  s^  ,  per  unit 

time,  including  the  heat  supply  arising  from  all  other  constituents,  with  r 

the  total  heat  supply  per  unit  mass  of  the  mixture;  the  internal  heat  supplies 

will  be  specified  by  constitutive  equations  but  do  not  affect  the  value  of 

r  in  (4.2),  .  Also  U  is  internal  energy  of  s  per  unit  mass,  allowing  for 
4  a  a 

interaction  between  s  and  the  other  constituents.  With  the  help  of  ( 3 -7 ) 
a 

the  point  form  of  (4.1)  becomes 

a 

v  D  U 

pr  -  q,  ,  -  Z  (o  -a  +  m  U  ) 
y  k,k  a-i  a  Dt  a  a 


_  ,  a  or  or  er  l  o  Of.  . 

+  E>i  Vi+°*i  Vi,k+2  maui  ui>  =  0  ‘ 

a  =  1  ’ 


(4.3) 


11. 


Except  when  all  the  temperatures  are  equal  we  cannot  develop  the 
consequence  of  the  inequality  (4.8)  with  the  help  only  of  (4.3).  We 
need  energy  equations  for  each  constituent  so  as  to  be  able  to  eliminate 
r^  frcm  (4.8)  and  also  to  provide  sufficient  equations  for  the  tempera¬ 
tures  T  (or  =  1 , . . , v)  .  We  consider  this  in  the  next  section. 
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5 .  Energy  balance  for  each  constituent 


Recalling  quantities  defined  in  previous  sections  we  postulate  an 
energy  balance  for  s^  in  the  form 


3  f  /lT  1  oi  or.  f  .  i  a  a.  a  f  1  A<*  Acr 

tt  p  (U  v.  v.  )dV  +  p  (U  tTV.  v.)v,  n,  aA-  tt  m  v.  v. 
3t  J,.  a  2  X  i'  JA  Har  a  2  i  i'  k  \  J,.  2  a  i  i 


dV 


r  ,  ,  O'  Of.  Of  Of  OfQf, 

=  ,,  [p  (r  +  F.  v.)  +  p.  v.  +  X,  .  ou.,j  dV 

JV  lKa  a  i  i  i  i  ki  ikJ 


J  <*1  Vi-  \  \)dA  +  L  *>a  dV  ■  JA  \  \  **  • 


(5.1) 


The  scalar  ^  represents  a  volume  contribution  to  the  energy  equation 

Of 

arising  from  interaction,  in  addition  to  interactions  already  considered, 

—of  —a 

and  n^  q^  represents  an  interaction  surface  flow  of  energy  with  q^ 

an  interaction  flux  vector.  However,  in  point  form  this  vector  contributes 

“Of 

q  which  can  be  absorbed  into  i)i  ,  so  without  loss  of  generality  we 

K  ,K  Of 

“Of 

set  q^  =  0  .  The  constituent  s^  has  a  mass  increase  acquired  from  the 
other  constituents  and  the  rate  of  change  of  its  kinetic  energy  just 
before  becoming  a  part  of  s^  is  represented  by  the  third  group  of 
terms  in  (5.1).  All  other  quantities  in  (5-l)  have  been  defined  pre¬ 
viously  and  follow  naturally  from  the  discussion  of  sections  3,4. 

Equation  (5.1)  can  be  reduced  to 
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O' 

D  U 


i  f  a  ..  /  f.ot  a.  a  1  a  a  1  Aa  Aa  , 

;  p  --  +  m  U  +  p(i.-F.)v.  +  -m  v.  v.  -  -  m  v.  v.  }  dV 
,,  a  Dt  a  a  Ha  i  i  i2aii2aiiJ 


_  r 


p  .  a  a  a  a  ,  p  ,  a  a  a> 

i  (caVPi  Vi+Xki  “ik+*„)dV  +  J„(ti  Vi_nk  qk)dA 


(5-2) 


1/,  1  AaAa  aa  a  a. 

I($+-^mv.  v.  +  p.  v.  rX.  .<d.,)  =  0  > 

a  2  a  i  i  ii  ki  ik'  ’ 


a=l 


(5-3) 


ana  if  we  also  use  (4.2),  then  by  summation  of  equations  (5.2)  we  find 
that 


or 

D  U 


[  Z  {p  7-a  +mU  -  p  r  +p(f?-F.  )va  +  ^  m  v?  v?}  dV 
JVa=i  a  W  a  a  a  1  1'  1  2  <*  i  i 


=  J  <  s.  h  -  \\)dA 

A  a=l 


(5.4) 


Equation  (5-4)  in  the  same  as  (4.1)  and  with  the  help  of  (3-7)  we  recover 
the  point  form  (4.3). 

Using  (3.8),  as  well  as  (3.13)  and  (3-3),  equation  (5*3)  reduces  to 


.aa  a  a  l  a  a  „  .  _ 

E  (p..  V.  +  Or  .  lU)..  +  -x  m  u.  u.  -  Y  )  =  0  , 

1  i  LkiJ  lk  2  a  i  i  a 


(5.5) 


where 


1  ,Aa  a,  .Aa  . 
Y  =  t  +  =■  m  (v.  -  v. )  (v.  -  v. ) 

a  a  2  j  i'  '  i 


(5.6) 


Also,  with  the  help  of  (3-7),  (3.13)  and  (5.6),  the  point  form  of  (5.2) 
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becomes 


P 


r 

at  a 


a 


D  U 

.  _ Of 

a  Dt 


m  U 

a  a 


«  rP 

a(ki)dik 


Y  =  0 

a 


(5.7) 


a  ot 

where  a^k^j  is  the  symmetric  part  of  cr^  .  The  arbitrary  functions 
’f  are  subject  to  the  condition  (5-5).  It  is  desirable  to  simplify 
this  condition  and  one  way  to  achieve  this  is  to  put+ 


Y 

a 


at.  a  v* 
a. (v.  -  v. } 
l  i  i 


a  .at 
+  CT[ki](u5ik 


v  .  i  a  a 
ui . ,  )  +  —  m  u.  u .  + 
ik'  2  a  i  i 


Pi 

a 


(5.8) 


so  that 


£  ©  =  0 
a=l  a 


(5-9) 


and  (5-7)  becomes 


Of  D  Ug  1  at  at 

pora  '  qk,k  *  pa  “DT  "  n«U«  +  2  ma  ui  ui  +  ®a 


at.  a  v.  a  a  at  .  at  v  .  . 

Hi(vi-vi)  +  a(ki)d.k  +  a[ki](U)ik  ‘  ®ik)  =  0 


(5-10) 


We  observe  that  for  some  purposes  it  is  useful  to  replace  the  v  equations 
(5.10)  by  an  equivalent  set,  namely  equation  (4.3)  and  the  v-1  equations 


+This  is  not  the  only  possible  form  for  Y  •  For  example,  the  first 

at  aa 

terms  in  (5.8)  could  be  replaced  by  (v“-v^)  »  the  difference  being 
absorbed  in  0  ,  but  (5-8)  is  more  convenient  for  our  purpose. 
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r  -  p  : 
a  a  v 


a  v 

qk,k  +  qk,k 


a  v 

DU  DU 

_ _ or  _ v 

pa  Dt  +  pv  Dt 


m  U  +  m  U 
a  a  v  v 


1  aal  vv  a,  a  v. 

+  m  u.  u.  -  77  m  u.  u.  +  0  -  0  +  u.(v.  -v.)  + 

S>  /kj  ^  n  i  i  Mi  .  .  “i't  ~i  ' 


d* 


:(ki)uik 


V  ,v  a  ,  a  V  . 

a(ki)dik +  ^kiiKk-^ik)  = 0  • 


(5.11) 


/v  a  a  ot 

In  equation  (5.10)  the  functions  5  qk  ’  ffla  are  aqq 

objective . 

If  we  eliminate  from  (4.8)  with  the  help  of  (5.10)  we  obtain 

the  inequality 


a 

D  A 


a 

D  T 


-  £  { p  (— -~  +  S  ■-)  +  mA  -  5-  m  u.  u.  -  0  } 

„  ,  T  lva  Dt  a  Dt  a  a  2  a  i  i  a 

a=l  a 


a  m 

v  T_ 


_  1  r  a,  a  vx  a  ,a  a  ,  a  v  xx  ^  a,k _ 

"aL^  ^  Vi_Vi  +^ki)  ik+<7[ki]  “’ik'<"ik  ’aflT'^  5  (5  } 


where 


A  =  U  -  T  S 

a  a  a  a 


(5.13) 


The  inequality  (5*12)  may  be  written  in  an  alternative  form  which  is 
useful  for  some  purposes .  Since 


17. 


a  a 

v  p  D  A  D  T 

a=l  a 

a 


m  A 

a  o 


a  Dt 


a 

A  D  T 


-  *  i°ah  &  *  4"  <s„  *  ?>  Tt2 ) 

a  a  a  a 


a=l 


da  j:  9  / 

“  P  Dt  +  E  3x  ^ 
0=1  k 


\ 


v  p  S 

)  *  I  -*■*  ( 

0=1  O 


DT 

_ g 

Dt 


3T 


’  <5-*> 


where 


S  =  S  +  ~ 
a  o  T 

o 


pA  =  E 


v  p  A 
o  o 


T 

o=l  a 


(5.15) 


the  inequality  (5.12)  becomes 


DA  *  P~S-  DT~  v 
P  nt  “  £ 

XJ  w  -.x  u  w  ■, 

0=1  O  0=1  o 


o  O  _  a  „  1  ,1  or  a  _  » 
T_  Dt  +_  T  %  ma  Ui  ui  +  ®a^ 


a  * 

v  p  u.  A  v  Q  T 

.  r  _a_  (  k  a  )  .  r  k  »>k 

o=l  To  ;  o=l  T2 

a 


*  1  r  O.  O  V.  0,0  0,0  V  ,  ,  , 

* *  -  W(vi-vi) +  ^(ki)dik +  a[ki)Kk-u,ik)}  - 0  ’  (5-16) 

o=l  a 


where 


*o  o  o  * 

^k  qk  +  po  To  ^  So 


(5-17) 


We  observe  that 


18 


o=l  o  a=l 


«'  k  '■  a,i  «  * 


Hence,  if  we  put 


90 

~  a  -a 


u“  =  T  t— ®  +  u,a  (a=  1, . .  ,v-l)  ,  a“-  =  T  0  6,  .  +  a*.  (o=  1, . .  ,v 

a  9x^  pi  '  ’  ’  '  ’  ki  r  a  ki  ki  ’ 

the  inequality  (5-16)  reduces  to 


DA 
-  P  Dt 


P  S 

o  or 

DT 

o 

T 

o 

Dt 

y,  o 

-(vi' 

vi)  + 

*o 

qk 

T 

a,k 

1  ,-ar,  a  v.  -or  , or  -or  ,  or  v 


-  Z 

a=l  T" 


^  ?0  , 


and 


v  o 


—a 


M-i  =  -  Z  ui  ,  Z  Or.  .1  =  0 

1  o=l  1  a=l  LK1J 


For  later  convenience  we  introduce  the  notation 


0 


o 


_1_ 

T 

o 


*o  = 


e  -  e 

a  v 


o  o  v 
Fik  "  “'ik’^ik  ’ 


zo  that  (5.20)  becomes 


(5-18) 


»  (5.19) 


(5.20) 


(5.21) 


(5.22) 


19- 


20. 


lonctitutive  equations 


We  assume  that  mass  elements  of  each  constituent  are  conserved 

o  that 


m  =  0 

a 


(6.1) 


ar.  s  that  the  constituents  of  the  mixture  have  elastic  and  viscous 

o  a  .  . 

nrorerties .  We  therefore  assume  that  A  ,  S  ,  o/,  q,  (or  =  1, . . ,  v) , 
•  -  a’  a  (hi)’  k  >  >  /> 

a  a 


,  .  o . - .  Z  (a =  1 , . . , v-l) ,  are  functions  of 
.  ki J  o'  ' 


j  re  5  xi 

axj  ’  ijk  sx$  axj 


(6.2) 


ana 


ep,k  ’  Vk 


VP 

i,k 


(6.3) 


for  g=l,..,v  ,  y=l,..,V“l  .  In  view  of  invariance  conditions  under 
superposed  rigid  body  motions  of  the  whole  mixture ,  the  velocity 
gradients  must  be  replaced  by 


(6.4) 


+ These  constitutive  equations  satisfy  equipresence  which  we  regard  as  a 
mathematical  convenience  and  not  a  physical  principle. 
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If  w e  use  the  inequality  (5.12)  to  obtain  restrictions  on  these  constitu¬ 
tive  assumptions,  it  can  be  shown  that  A  depends  only  on 


T 


0 


(6.5) 


In  addition,  further  restrictions  are  placed  on  A  and  other  dependent 

a 

functions,  but  there  is  some  algebraic  complexity  in  making  explicit 

deductions  about  the  form  of  A  from  these  restrictions. 

a 

Inspection  of  (5.20)  suggests  an  alternative  procedure.  We  assume 

the  Aa  ’  Sa  (and  henCC  A’  S’  0ty)’  ®(ki)’  %  4®  > 

a“k.],  0^  (a=  1, . .  ,v-l)  are  functions  of  the  variables  (o.2)  and  (6.4). 

Prom  the  inequality  (5.20)  we  again  find  that  A  depends  only  on  the 

variables  (6.5).  From  (5.17),  (5*l8)  and  (5.19)  we  see  that  q*  and 

aki  are  functions  of  the  variables  (6.2)  and  (6.4),  but  that  pk  are 

functions  of 


%kr 


9G? 

ilk 


ax; 


P 


ik,r 


r3 

1  ik,r 


as  well  as  the  variables  (6.2)  and  (6.4).  If  we  add  the  additional 

restriction  that  p.®  should  depend  only  on  the  quantities  (6.2)  and 

(6.4),  it  follows  that  0  reduce  to  functions  of  the  variables 

a 


9 

v 


("•6) 


Using  (5.18)2  and  the  fact  that  A  is  a  function  of  the  quantities 
(6.5),  it  follows  that  A^  reduce  to  functions  of  the  variables  (6.6). 
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■vistri  -tions  on  the  remaining  constitutive  equations  may  now  be  found 
from  (5.12),  (5.I0),  (5.20)  or  (5.23),  and  we  use  (5.23).  This  procedure 
appears  to  be  slightly  more  restrictive  than  the  first  and  is  the  one 
adopted  here. 

We  aaopt  the  notation 

h  ■  tv  Vk-  dTk'  O  (A.  1, . .  ,l6v-7)  ,  (6.7) 

and  we  put 
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qk 

0 qk 
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*Or 

*a 

qk 
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eqk 
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„0 
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a 

0  or 

e  a 

(6.8) 


where  the  quantities  with  a  prefix  o  are  the  values  of  the  corresponding 
quantities  on  the  left  hand  side  of  (6.8)  when  Y^  =  0  .  The  remaining 
terms  in  (6.8)  may  be  regarded  as  polynomials  of  degree  greater  than  zero 
in  the  functions  Y^  with  coefficients  which  are  functions  of  the  remaining 
variables  in  (6.2). 
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With  the  help  of  (5-20)  we  have 


*  _  a  BA  _  ~  BA- 

pcrSa  '  p9a  30  pTo  3T 

a  ot 


(6.9) 


so  that 


v  BA 

p  s  =  e  e  paefl  =  - 

«  «  «  p»i  6  *  aea 
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tK  „  ,  T„  5T 
0=1  0  a 


(6.10) 


Also,  we  write 


A'  +  A" 
a  a 


A  -  A'  +  A"  ,  pA'  =  6  E  p  A' 


v«.l  * 


(6.11) 


where  A1  are  the  values  of  A  when  A-  ,  V?  vanish,  and  A  is  the 

a  or  P  k 

value  of  A  when  Ag  vanishes.  For  the  present  we  regard  a"  as  a 
polynomial  of  degree  greater  than  zero  in  Ag  ,  and  a"  as  a 

polynomial  of  degree  greater  than  zero  in  Ag  with  coefficients  which 
are  functions  of  6  ,  F*  .  The  result  (6.11k  follows  from  (5.15)2- 

V  1 J 

Again,  using  (5.23)>  we  have 


-a  _  _£_  BA'  F<* 

ocki  "  0  kj 

v  BF.^ 


(6.12) 


-Of 

o^k 


1  '  .  _BAl  _0  n  .^Al  G«  ^jr  ) 

■  e.  a^(p«  iJr  axP  Pe  iJ1-  >v° 


v  0=1 

for  or  =  1, . .  ,v  provided 


3FiJ 


(6.13) 
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(6.14) 


Jl  3f“.  * 

ij 


3Fkj 


Also 


P 

o  a 


0 


0 


(6.15) 


xn  audition  there  is  a  residual  inequality  which  we  omit.  From  (5.18), 
(5.19),  (6.12)  and  (6.13)  we  then  have 
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o°ki 
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=  I 
3=1 


p3  dr  Fkj 
aFij 
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{i,  -  I!  (p 
°  k  P=1  *  3F7. 
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3A '  .  3XP 

-  g  r3  _ £  .  „  _ w  r- 

3  ijr  P  p3  or  Gijr  or 

ij  3  k  aFij  Sxk 


(6.16) 


(6.17) 


in  agreement  with  previous  results  [9] . 

The  above  constitutive  equations  are  subject  to  the  usual  invariance 
conditions  under  superposed  rigid  body  motions  of  the  whole  mixture. 

For  example,  A  ,  A  reduce  to  the  new  forms 


A 
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3xa  BxP 


u**  U-'*.  dx  Q 

A  (e  ,  a  ,  — E  ,  — E  v*  ) 


3X7  3X? 
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a  g 
3x  9xp 

A  =  A(0  ,  A  ,  — -  — -  ) 
v  v  a  e  ' 
Y  3X.  3X . 

^  J 


(6.18) 

(6.19) 


for  Qr,0=l,..,v  ,  \  =  1, . .  ,v-l  .  in  view  of  (6 .19)  we  may  verify  that 
the  condition  (6.14)  is  satisfied. 

In  obtaining  the  above  constitutive  equations,  we  have  adopted  the 
second  procedure  described  earlier  in  this  section  [after  (6.5)].  We 


25. 


note  that  if  we  had  used  the  first  method,  we  would  still  recover  the 
results  (6.9),  (6.12),  (6.13),  (6.l6)  and  (6.17)  with  being  the 
"equilibrium"  values  of  A^,  i.e.,  the  values  of  A^  when  vanish. 

The  special  case  of  the  formulae  (6.10)  when  all  the  temperatures 
are  equal  is  of  seme  interest.  In  order  to  see  what  happens  to  (6.10) 
we  write 


0  =  A 
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0  -  0  =  A  (0=  1, . .  ,\j-1) 
a  v  a 


then 


(6.20) 


and 


IT -ST  • 

a  a 

_3 _ 3_  V‘1  _b_ 

90  ~  9A  9A 
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Qf=l 


V-l 

+  Z 
a=»l 


or  a  a. 


(6.21) 


(6.22) 


The  functions 


may  be  expressed  in  the  form 


V  Ae(Fij’  ^  V  •  (6>23) 

When 

Ay  =  0  (Y-l,..,v-l)  ,  A  =  1/T  ,  (6.24) 


where  T  is  the  common  temperature  of  the  constituents ,  we  assume  that 
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*.  a 

VFij’ 


(6.25) 
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(6.26) 


BA 


BA 


.2  BT 


It  will  not,  in  general,  be  possible  to  evaluate  BA^/BA^  as  A  -•  0 

* 

in  terms  of  the  function  A  and  we  assume  that  these  derivatives  are 

P 

replaced  by  arbitrary  finite  scalar  functions^  From  (6.10)  and  (6.21) 
we  see  that  under  the  condition  (6.24), 
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v  v  o  BAR 

I!  p  S  =  T.  A  pfl  vr~ 
a=l  “  a  p-1  * 


-  E  P 
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(6.27) 


or,  defining 
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pS  =  £ 
Qf=l 


P  S 
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pA  =  S 

P-1 


V,* 


(6.28) 


then 


S  =  - 


BA 

BT 


(6.29) 


as  given  previously  [6,8] .  Also  S  -S  become  arbitrary  functions  of 

a  v 

6  B  * 

F. . ,  V,  ,  T  not  expressible  in  terms  of  A_  ,  which  again  agrees  with  a 
K  P 

previous  result  [9J.  We  recall  that  the  energy  equations,  which  provide 
v  differential  equations  for  the  v  temperatures  ,  may  be  taken  to 
be  (4.3)  and  the  v-1  equations  (5.11),  and  we  observe  that  volume  supplies 
of  heat  may  be  taken  to  be  given  by  the  v  independent  quantities 


+This  is  analagous  to  the  procedure  used  in  the  passage  from  a  compressible 
to  an  incompressible  elastic  material. 
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(e=  i,..,vi) 


(6.30) 
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pr  =  E 
or=l 


P  r 
or  o 


pere  ' 


p  r 

V  V 


When  temperatures  reduce  to  a  single  temperature  we  have  one  equation 
(4.3)  for  this.  The  remaining  equations  (5.11)  can  then,  in  general, 
only  be  satisfied  by  appropriate  choices  for  the  heat  supplies 
p.r  -  p  r  .  This  is  similar  to  the  situation  in  which  we  can,  in 

p  p  V  V 

general,  only  maintain  isothermal  conditions  with  T  constant  by  an 
appropriate  choice  of  r  . 

Special  cases  of  the  results  (6.16)  and  (6.17)  may  be  obtained 

without  difficulty.  For  example,  if  constitutive  equations  depend  on 

the  deformation  gradients  F? .  only  through  the  densities  p  ,  then 

ij  p 


a  "_B 

oCTki  "  "  6ki  papB  Bpa  ’ 

of  *  ,  dpa  v 


(6.31) 


(6.32) 


as  found  previously  [6],  These  results  are  appropriate  for  a  mixture  of 
viscous  fluids. 
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